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Abstract We present a global study of low mass, young stellar object (YSO) surface 
densities (E) in nearby (< 500 pc) star forming regions based on a comprehensive 
collection of Spitzer Space Telescope surveys. We show that the distribution of YSO 
surface densities is a smooth distribution, being adequately described by a lognor- 
mal function from a few to 10^ YSOs per pc^, with a peak at ^ 22 stars pc^^. The 
observed lognormal 1, is consistent with predictions of hierarchically structured star- 
formation at scales below 10 pc, arising from the molecular cloud structures. We do 
not find evidence for multiple discrete modes of star-formation (e.g. clustered and 
distributed). Comparing the observed Z distribution to previous I, threshold def- 
initions of clusters show that they are arbitrary. We find that only a low fraction 
(< 26%) of stars are formed in dense environments where their formation/evolution 
(along with their circumstellar disks and/or planets) may be affected by the close 
proximity of their low-mass neighbours. 



1 Introduction 

It is often stated that most if not all stars form in clusters. This view is based largely 
on near-infrared (NIR) studies of star-forming (SF) regions within several hun- 
dred parsecs of the Sun ll27]|29l . However, combining high-resolution mid-infrared 
(MIR) data with the NIR makes YSO identification more robust and less likely to 
suffer from field star contamination, which leads to better tracing of YSO surface 
densities. This means that with the NIR alone, there were large uncertainties in the 
number of stars at low values of YSO surface densities (Zyso) El- 

With the launch of the Spitzer Space Telescope [351 we are now able to differen- 
tiate YSOs and contaminating sources based on colour information and hence can 
study the distribution of YSOs independently of the surface densities. Large field- 
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of-view (FoV) Spitzer observations of SF regions ||T] [16] found that YSOs extend 
well beyond the densest groups in their environment and continue throughout. Sev- 
eral Spitzer surveys that cover nearly all the SF regions within 500 pc of the Sun 
were combined and presented by 17] . 

The spatial distribution of forming stars, i.e. do they form in clusters, is important 
for many reasons and here we list several. What spatial distributions of star forming 
environments lead the Galactic stellar clusters and why do they constitute such a 
low fraction of the stellar population? Are there multiple discrete modes of star 
formation, such as clustered and distributed, that manifest themselves as peaks in 
a surface density distribution (e.g. Il32l [8] |34] |33])? What fraction of star forming 
environments are dense enough to affect the YSOs and alter their disk and planet 
formation/evolution |[T]|20l, which we refer to as dense environments. 

Outside the Milky Way we see that star-formation proceeds hierarchically from ^ 
10 pc (the resolution limit of most extragalactic surveys) to kpc scales (e.g. ITOirTTl 
mis] [it])- If star formation occurs hierarchically down to 0.1 pc, then it is unlikely, 
and indeed not expected, that all stars form in gravitationally bound young stellar 
clusters. The availability of large FoV Spitzer surveys of nearby star-forming regions 
allows us to test how far down in scale this hierarchy proceeds. By hierarchical 
structure we mean a smoothly varying non-uniform distribution of densities, where 
denser subareas are nested within larger, less dense areas [BTlfT?! . 

In this paper we will investigate 1 ) whether there is evidence for multi-modality 
in the surface densities of YSOs, 2) what fraction of stars form in dense environ- 
ments in the local neighbourhood and 3) how relevant the various cluster definitions 
are as provided by i] |27] [D ED E] . 



2 Data & Biases 

Multiple Spitzer surveys were used to generate a comprehensive and statistically 
significant dataset to investigate the spatial surface density properties of forming 
stars in the solar neighbourhood. The surveys are the Gould's Belt (GB) survey 
(Allen et al. in prep.), Orion survey (Megeath et al. in prep.). Cores to Disks (c2d) 
survey ifTSI . and the Taurus survey ll30l . The GB and Orion catalogs have not been 
publicly released yet. There are more than 7000 YSO detections in the combined 
catalogs at distances between 100 to 500 pc. 

Spitzer data are necessary for this study as low i^yso can be differentiated from 
field star populations, unlike NIR observations where field star contamination can 
be problematic. The comprehensive YSO population represents a global view of the 
low-mass star-forming region in the local neighbourhood from low to high surface 
densities. These Spitzer surveys combined represent the most complete census of 
star formation within 500 pc of the Sun available to date. 

Data treatment on distances, exclusion of the central region of the Orion Nebular 
Cluster, rejection of Class 111 objects, and potential contaminants are described in 
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Fig. 1 (a) The cumulative fraction of surface densities for the GB+Taurus, c2d, and Orion sur- 
veys. Each SF region included in the distributions has N(YSOs) > 10 and a sufficient field-of-view 
to properly calculate stellar surface densities. The Orion survey stops at 73% for the cumulative 
fraction since the ONC is excluded. We adopt a 65% disk fraction for all of the SF regions. We 
normalised each curve by the number of YSOs in each survey, (b) With the GB+Taurus, c2d, and 
Orion surveys combined we see Class I & II distributions having similar profiles with a small off- 
set in density, showing that we are likely seeing the primordial distribution of the YSOs. (c) With 
all of the Spitzer surveys combined we compare several cluster definitions. The vertical grey lines 
from left to right are Lada & Lada (2003), Megeath et al. (in prep.), J0rgensen et al. (2008), Car- 
penter (2000), and Gutemiuth et al. (2009) stellar density requirements for clusters. These values 
correspond to 3, 10, 20, 32, and 60 YSOs pc^^ and intersect the corrected cumulative distribution 
profile, implying that 87%, 73%, 62%, 55%, and 43% of stars form in clusters, respectively. The 
percentages coirelate to what fraction of stars foim in 'clusters' based on the various definitions. 
The black vertical line is for a dense environment where Z > 200 YSOs/pc^. The fraction of YSOs 
in a dense environment is < 26%. 
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Q. After the data correction process, we have 3857 YSOs remaining to estimate 
the surface densities presented in this paper. 



3 i^YSO distributions 

Our primary tool for analysing the surface densities is computing the local ob- 
served surface density of YSOs centred on each YSO's position, where Eyso = 
(N — l)/(7rD^) and N is the Nth nearest neighbour, and Dn is the projected dis- 
tance to that neighbour (see |[9l)- Throughout this work we will adopt N = 7, al- 
though we note that all results have been tested for N =4 — 22 and no significant 
differences were found. Figure [T]shows the surface density distribution of all YSOs 
in our sample. 

In order to see the fraction of YSOs above a given E threshold, we show the com- 
bined i^YSO distribution (shown as a cumulative fraction normalised to the number 
in each combined survey) for the three surveys used in this study in Fig. [T^. Note 
that the GB/Taurus distribution lies to the left of the c2d survey. This is simply due 
to the GBATaurus focussing on lower density regions than c2d. The cumulative dis- 
tribution for the Orion survey only reaches 0.73 in Fig. [T^ and 0.8 1 in Fig.[T];, where 
all the surveys have been combined, due to the exclusion of the ONC. In Fig.[T]; we 
show the cumulative distribution of all YSOs included in our survey, while in Fig.[TJ) 
we split the survey into class I and class II objects. 



4 Results 

It has been long assumed that two distinct modes of star formation exist for YSOs, 
'clustered' and 'distributed' (e.g. |fT9l H] EtJ), but the notion has been questioned 
after Spitzer results hinted otherwise |[T]. If there are indeed two modes, then we 
would expect to see a bi- or multi-modal profile in cumulative surface density dis- 
tribution plots such as Figs. [Th ands & [TJ;. Instead we see smooth and featureless 
distributions from the low to high stellar surface densities for the c2d, GB, Taurus, 
and Orion surveys. We find that the Zyso distribution of low-mass stars in the solar 
neighbourhood can be well described by a lognormal function, as seen in|2] with a 
peak at ^ 22 YSOs/pc^ and a dispersion OiogigX = 0.85. 

The spatial distribution of the YSOs in these SF regions is expected to be close to 
primordial since their YSOs, in particular Class Is and Clas s lis, are < 2 Myr old 1221 
l24ll . In order to place stricter constraints on this, we now split the complete sample 
into Class I and II objects, which can be roughly attributed to an age sequence. 
The cumulative E distributions of Class I and II YSOs are shown in Fig. [TJ). We 
see that the two distributions have similar smooth density spectra, however they are 
slightly offset. The E of the Class I/II objects are calculated by finding a YSO's Nth 
nearest YSO. Once this is done for the YSOs we separate the Class I/II objects. E 
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Fig. 2 Top panel: The surface density distribution of the total sample of YSOs in the solar neigh- 
bourhood used in this work (black). A lognormal function with a peak at ~ 22 YSOs/pc^ and a 
dispersion ffiogigZ = 0.85 is shown as a dashed line. Bottom panel: The same as the top panel but 
now broken into the three respective surveys. 



is calculated this way since Class Is and Class lis are not always spatially distinct 
from one another ED . Class lis are known to be slightly more dispersed than Class 
Is in high density regions 1211 reflecting early dynamical evolution. However, the 
similar distribution between these classes leads us to conclude that the distribution 
of observed L is mainly primordial in nature. 

In Fig.[l]: we show five vertical grey lines that refer to the defined densities re- 
quired for a collection of YSOs to be considered 'clustered' (Lada & Lada 2003; 
Megeath et al. in prep.; J0rgensen et al. 2008; Carpenter 2000, Gutermuth et 
al. 2009). The vertical lines fall on the same featureless slope and do not corre- 
spond to any preferred density. The black vertical line, which corresponds to dense 
environments (as defined in ll20l ). shows that < 26% of YSOs are formed in envi- 
ronments where they (along with their disks and planets) are likely to interact with 
their neighbours. 



5 Discussion and Conclusions 



We conclude that stars form in a broad and smooth spectrum of surface densities 
and do not find evidence for discrete modes of star formation in the E of low mass 
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YSOs forming in the solar neighbourhood. The observed lognormal surface density 
distribution is consistent with predictions of hierarchically structured star-formation, 
which arise from similar structures in molecular clouds. ( [fT2l[T4l ). Our results sug- 
gest that clusters are not fundamental units in the star formation process, but simply 
the high density tail end of a continuous distribution. A small fraction (< 26%) of 
stars form in dense environments where their formation and/or evolution is expected 
to be influenced by their surroundings. 

By comparing our global surface density distribution with the summary of clus- 
ter definitions in (71, we find that the fraction of stars forming in clusters is crucially 
dependent on the adopted definitions (ranging from ^ 40 to 90%). Since star/cluster 
formation happens in a dynamic environment ^ |2] [28] it may be impossible to 
define what will become a stellar cluster in these early stages, as the final object 
that remains bound is simply the dynamically mixed part of a larger, initially hi- 
erarchical, distribution. In this scenario a cluster can only be clearly defined above 
the surrounding distribution once it is dynamically evolved where tage/tcross > 1, as 
defined by Gieles & Portegies Zwart 2010. 

Star forming environments provide the initial conditions from which star clus- 
ters may eventually form, albeit rarely. Since the probability density function of 
molecular gas varies with environment, as does the tidal field experienced by the SF 
regions, it is likely that the fraction of YSOs ending up in bound star clusters varies 
with environment lfT4l and the observed Eyso is not universal. Future investigations 
will extend this work out to 2 kpc where proto-stellar cores and high-mass SF re- 
gions will be included, and testing whether star formation primarily occurs when 
Zgas is above a specific threshold (e.g. ll23l [3] |26l ) 
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